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Abstract 

The Ising approximation of the Hcisenberg model in a strong magnetic field, with two, three and six spin 
exchange interactions is studied on a kagome chain. The kagome chain can be considered as an approximation 
of the third layer of ^He absorbed on the surface of graphite (kagome lattice). By using dynamical approach 
we have found one and multi-dimensional mappings (recursion relations) for the partition function. The 
magnetization diagrams are plotted and they show that the kagome chain is separating into four sublattices 
with different magnetizations. Magnetization curves of two sublattices exhibit plateaus at zero and 2/3 of 
the saturation field. The maximal Lyapunov exponent for multi-dimensional mapping is considered and it 
is shown that near the magnetization plateaus the maximal Lyapunov exponent also exhibits plateaus. 
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1. Introduction 

The partition function of a statistical system can be calculated by means of several methods: transfer 
matrix approach dynamic technique (H-Q etc. The most useful for recursive lattices is the dynamical 
approach, which is a very powerful tool for research of a number of theoretical problems in statistical 
mechanics. By means of this approach one can obtain exact recursive relations for the partition function on 
recurrent lattices. 

As the systems of almost localized identical fermions one can model solid and fluid ^He films. Since 
the light mass spin-1/2 ^He atoms are subject to a week attractive potential, the theoretical explanation 
of magnetism is based on the multiple-spin exchange mechanism (5|. An important case is represented by 
solid and fluid ^He films absorbed on the surface of graphite 0, since it is a typical example of a two- 
dimensional frustrated quantum-spin system Q . The first and second layers of the system form a triangular 
lattice, while the third one forms a system of quantum 1/2 spins on a kagome lattice [lo| . 

In a strong external magnetic field Heisenberg model can be approximated by the Ising one. If the strong 
magnetic field is directed along the z-axis, a reduction of transverse fluctuations is expected to occur; it 
is supposed that in this case and cr^-spin components are infinitely small infinitely small and may be 
neglected 

The phenomenon of magnetization plateau has been studied during the past decade both experimentally 
and theoretically. The plateaus may be exhibited in the magnetization curves of quantum spin systems at 
very low temperatures in case of strong external fields. The plateaus of this kind were first observed by Hida 
[l^ in trimerized S=l/2 spin chains and they were analytically explained by Okamoto [l3| . Magnetization 
plateaus appear in a wide range of models on chains, ladders, hierarchical lattices, theoretically analysed 
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by dynamical, transfer matrix approaches and exact diagonalization in clusters (see Ref. |14|-j2a|). The 
experimental evidence for the appearance of magnetization plateaus has also been found for two-dimensional 
systems l2 7| - [ 3Cl . On triangular lattices a magnetization plateau was observed at m/ms = 1/3 for compounds 




^sCuCls 
H)- 



[32| (see also Ref. [33[-[35|); plateaus have also been observed on a kagome lattice 



like CgEu 
(see Ref. 

In the present paper the dynamical (recursive) approach was used to study magnetic properties of a 
kagome chain. The magnetization plateaus were found at low temperatures and moreover, the kagome chain 
was separated into four sublattices with different magnetizations. Two of these exhibited plateaus, whereas 
the others did not. The stability of the system is described by Lyapunov exponents [41 by means of which 
the exponential rate is measured, at which the adjacent orbits converge or diverge |42| - 4J]. It is interesting 
to check whether the maximum Lyapunov exponent has plateaus that coincide with the magnetization one. 

This paper is organized as follows: in Section 2 the two, three and six-site exchange interactions in a 
strong magnetic field are described, in Section 3 the exact recursion relation for the partition function is 
derived, the magnetization curves for different values of temperature and exchange parameters arc displayed 
in Section 4. In Section 5 a multidimensional mapping for the partition function is obtained and the 
corresponding maximum Lyapunov exponent is investigated. Finally, Section 6 contains the concluding 
remarks. 



2. Exchange Hamiltonian wiih two, three and six- site exchange 

The Heisenberg Hamiltonian for kagome chain consists of two parts 

H = H,,+Hz, (1) 

where H^x is the spin exchange Hamiltonian and Hz is the Zeeman Hamiltonian which is responsible for 
magnetism, which has the following form 

Hz = -J2 i^^*^' = E (2) 

i i 

where 7 is the gyromagnetic ratio for a nucleus, and B is the magnetic field 
According to Q the multiple spin exchanges Hamiltonian can be written as 

He. = J2 ^ - J3 E (^3 + ^3"') +J6 E (^6 + Pe') , (3) 

Pairs Triangles Hexagons 

where P2 is the pair transposition operator, P3 (Pg) is the operator making a cyclic rearrangement in the 
triangle (hexagon). The explicit expression of pair transposition operator was been given by Dirac 

^2^P,, = i(l + T,^,), (4) 

where cr^ is the Pauli matrix, acting on the spin at the i-th site. The expressions for P3, and Pg in terms of 
Pauli matrixes were derived in (23j . By using these expressions one can rewrite the Heisenberg Hamiltonian 
in terms of Pauli matrices. By taking into account that in a strong external magnetic field directed along the 
z-axis the contribution from x and y spin component will be negligible and the only relevant contribution 
will come from the z component (which can effectively take values Sz = il), we may consider an Ising model 
instead of Heisenberg one. By using expressions for P3, and Pg and by substituting the Pauli matrix with 
Sz we get the Ising approximation of the Heisenberg Hamiltonian: 

^ = y E (i + ^'^^)-y E i^ + ^^^, + ^j^k + s,s,) + ^-^ (5) 

Pairs Triangles Hexagons i 
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Figure 1: Sublattices on the kagome chain. 

where represents the six particle exchange term in each hexagon: 

^(6) 



y li<v f2<iy<\<p 



SnSi^S^Sp + S1S2S5S6S9S10 



(6) 



where the first sum goes over all pairs in the hexagon, and the second one goes over all quartets in the 
hexagon. 



3. The recursion relation for the partition function 

The partition function of the system is defined in the following way 

Z = g-ff(si,S2,...S„)/fcBT 



(7) 



where H (si, S2, ■■■Sn) is the Hamiltonian of the system, ks is the Boltzmann constant, T is the temperature 
and the sum goes over all {si, S2, •■•Sn} configurations of the system. Below all the parameters will be 
rescaled using the Boltzmann constant (J2/fc_B = </2, Ja/^s = J31 h/fc^ = h). To obtain recursion relations 
for the partition function one can separate the kagome ladder into two identical parts (branches) and firstly 
perform a summation over all spin configurations on each branch, and secondly sum over the central spin 
variable (see Fig. 1). The summation on each branch yields the same result and such a term only depends 
on the value of central spin: 



(8) 



where gn{so) denotes the contribution of each branch. The expression for _g„(so) can be written (si . . . sn, 
see Fig. 1) in the form: 

S„(so)- E e'=(^--^")g„_i(sn), (9) 

Sl,S3---Sll 

where and n is the number of generations and the expression for k(so, ...sn) has the following form: 



fc(so,-Sii) = {l + SiSj) + ^ Y {i + s,Sj + SjSk + SkS.,) ~ ^ ■ H^^'> 

18 pairs 6 triangles 

Since gn-i does not depend on sqi ■•■SiO: equation Q may be obtained as 



T 



E*- (10) 



gn{+) = /v(+,+)g„-i(+)+X(+,-).g„-i(-), 
5„(-) = /v(-,+)g„_i(+)+X(-,-)g„_i(-), 



(11) 




-3-2-10123 -3-2-10123 



a) h b) h 

Figure 2: 

Magnetization curves for J2 = 2mK, J3 = 1.5mK, Jq = 2mK, T = Q.QAmK for a) different sublattices 

b) the average magnetization of all sublattices 



where we denote gn{si) by (7n(±) according to the sign of Si and if(so,sii) is the sum of fc(so,...sii) over 
the ten intermediate spins between si and S12 

X(si,si2)= e^(^«'-^"). (12) 

The recursion relation is obtained by introducing a new variable x„ = ): 



" i^(-,+).T„_l+i^(-,-)- 



(13) 



4. Magnetization plateaus 



The Thermodynamic functions of the system, such as magnetization, can be expressed in terms of Xn- 
The local magnetization at site is expressed as 

m^{s^) = ^. (14) 

where 

P{sr) = ^ s,e-^(^i'*^--^")/'=«^. (15) 
{^.} 

By using equation (|5]) one can write for the central vertex 

e*5^(+) + e-T52(_) eTx2 + 1 

To obtain local magnetization for non central vertices one should replace one of two multipliers (?n(si) in 
dHl) by formula ® (see Fig. 1) 

^ e*^° ■e'=(^--^") •5n(so)-5„-i(sii). (17) 

So,S2---Sll 
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Suppose it is required to find the local magnetization for the r-th site (?' 7^ 0, 11): in that case we must 
expand sums over sq, sh and Sr 

Z = (At,+ + !<+, + ) ■ 9ni+)gn-l{ + ) + + + KZ^+ + Kl_ + X- _) • 9ni+)gn-l{-) + 

+ {K+^_ + KZ,S)-gn{-)9n-i{-). (18) 

the sum 

Kll,^,= e*^«-e^(^--") (19) 

so..s,,_i,sv+i...sii 

is denoted as _(_ depending on the sign of sq: sh and sign. In the same way the expression for P{sr) 
may be obtained. The difference is that each ^ is replaced by {~K^ _|_): 

P{sr) = {kX,+ - K+^+) ■ 5n(+).9«-i(+) + {K-,+ - + /^|,_ - A7__) ■ 5„(+)r7„_i(-) + 

+ (A'+ _-/!-_) •5n(-).9„-i(-). (20) 

By using equations (|18p and (j20p we find the following expression for the magnetization at r-th site in terms 
of Xn 

~ (A-+ + + +) • XnXn-1 + (A'+ + + KZ^+ + A'+ _ + K- _) ■ x„ + {K+ _ + /fl _) ' 

For fixed values of temperature and exchange parameters the dependence of magnetization on the external 
magnetic field may be obtained by means of a simple iteration scheme induced by the recursion relations, 
starting from some initial value xq, The thermodynamic limit then corresponds to the asymptotic limit 
{n — !■ 00). 

Because of the symmetry of the Hamiltonian, the values of magnetization for some spins are equal: we 
have in general (See Fig. 1), 

771(53) = 774(54) = ^77(57) = 777(S8) — A 

m{si) = m{s2) = m{sg) = 7ri(sio) - B 
777(55) = 777(55) - c 

777(so) = 777,(511) - D (22) 

In Figure 2 magnetization curves are plotted for different sites. As it can be seen from the figure, the 
magnetization functions have plateaus only for the sublattices B and C. 

5. Two dimensional mapping and Lyapunov exponent 

It is interesting to calculate the Lyapunov exponent near plateaus but one cannot use recurrent function 
(Uni) because it is related to 5i site. To obtain the recursion relation for the partition function that is related 
to 52 site one can separate the kagome chain into two identical parts as shown in Figure 3: 

^= E ^A\p,^-<?n(si,S2)-5„(pi,P2), (23) 

where 

^Si,S2 ^ y^g{T(-'Sl+S2+^0+Pl+P2) + ^^^^(g<^''"(^0,.'il..'i2)+g'^'^'(so,pi,P2))-#g"''} (^2A) 
So 

and the contribution of each branch is denoted by g„(5i, S2). The recursion relations for can be obtained 
by performing eleven (53 . . . 513, see Fig. 3) summations in the kagome chain: 

3n(si,S2)- ^Ii'2,'?3 ■5n-l(si2,Sl3), (25) 

S12-S13 
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Figure 3: The kagome chain. 
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where K^}' is as follows 
* 12 5*13 



Si2,Sl3 



.73 - J2 

6 Triangle z = 3 



52 J£^(6)_^^ ^ 



(26) 



Since i^lj^^^l^^ = K^^ H^, K^^^ll^ = K^^^^/^^ one can show that .g„(+,-) = ,g«(-,+) and, hance, the recur- 
rence relations are two-dimensional. By introducing a;„ = .g„(+, +)/gn{—, ~) and ?/„ = ^)/5n(~i ~) 
the recursion relations may be written in the following form: 



x„ = $(a;„_i,2/„_i), ^ix,y) 
?/„ = ^'(x„_i,y„_i), 



A"+;+a; H- 2K~'~y + KZ'~ 



The magnetization can be expressed in terms of a;„,y„ with due regard for (|23p and (|24p 
m(si) 



E Sl ■^pi,^ ■5n(si,S2) ■ 9n{Pl,P2) 

Sl ,S2,pi ,P2 



Sl,S2,Pl,P2 

Taking into account that N^^'^^ = Np^'^l = iVp,^pJ = NP^f^, we get after the expansion of sums: 

2N+--yn - NI-Z 



m{si) 



Nt:txl + 4.N+'+Xnyn + iN+'-yl + 2N+'+Xn + 47V+Cy„ + NZ:Z 



(27) 



(28) 



(29) 



Magnetization curves calculated by formulae ((29)) and ([2T|) with r = 1 coincide. Relations ([27)) permit the 
calculation of Lyapunov exponents near the plateaus: in the two-dimensional space they have the form 



Ai,2 = lim ^\n{eigenvalues{A{xi,yi) ■ A{x2,y2) ■ ■ ■ MxN,yN)}) , 

N^oc iV 

where K{xi,yi) is Jacobian matrix evaluated at the [xi^yi) point 

/9* 



dy , 



(30) 



(31) 



The only subtle point in the estimation of Lyapunov exponents j45l |46| . |47| is in the rescaling procedure 
in each iteration to avoid the numerical overflow. The results of calculation are shown in Fig. 4 with the 
corresponding magnetization function: it is seen in the figures that the maximum Lyapunov exponent also 
exhibits plateaus. Moreover, for the maximum Lyapunov exponent the location of magnetization plateaus 
coincide with those of magnetization curves. 
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Figure 4: Coincidence of magnetization plateaus and the maximum of Lyapunov exponent. 



Conclusions 

In the present paper the theory of dynamical systems has been used to study magnetization on the kagome 
chain. In a strong external magnetic field the Ising model has been considered instead of the Heisenberg one 
for the case of two, three and six-site exchange interactions. Exact recursion relations have been derived for 
the partition function. The magnetization curves for different temperatures and exchange parameters were 
obtained and the kagome chain was observed to split into four sublattices with different magnetizations. 
Two of them exhibit plateaus at zero and 2/3 of the saturation field value at low temperatures and certain 
exchange parameters. A multidimensional mapping was deduced to study the maximum Lyapunov exponent 
near the plateaus and it was shown that the location of plateaus on the maximal Lyapunov exponent and 
for magnetization was the same. 
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